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Abstract 

The kinetic energy of bulk relativistic plasma ejected from the central engine of 
blazars is converted into nonthermal particle energy in the comoving frame through 
a process of sweeping up material from the surrounding medium. The resulting de- 
celeration of the bulk plasma introduces a number of effects which must be included 
in blazar modeling. For example, the varying Doppler factor means that model fits 
must employ time integrations appropriate to the observing times of the detec- 
tors. We find that the ratio of the peak synchrotron fluxes reached at two different 
photon energies is largest when viewing along the jet axis, and becomes smaller 
at large angles to the jet axis. This effect is important in studies of the statistics 
of jet sources. Variability due either to bulk plasma deceleration or radiative cool- 
ing must be distinguished in order to apply recently proposed tests for beaming 
from correlated X-ray and TeV observations. The blast-wave physics developed to 
analyze these problems implies that most of the energy injected in the comoving 
frame is originally in the form of nonthermal hadrons. Because plasmoid deceler- 
ation can produce rapid variability due to a changing Doppler factor, arguments 
against hadronic blazar models related to the long radiative cooling time scale of 
hadrons are not compelling. 
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1 Introduction 

Ejection of relativistic plasma from a compact central engine is thought to 
account for the appearance and observational properties of a number of fasci- 
nating systems in astronomy, including galactic black hole jet sources, radio 
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galaxies, quasars, and gamma-ray bursts (GRBs). Several arguments have led 
to this conclusion, perhaps the most important being measurements of appar- 
ent transverse superluminal motion in multi-epoch VLBI observations of radio 
quasars at the sub-pc scale. Apparent transverse speeds exceeding ~ 10c are 
found in many sources (e.g., Vermeulen & Cohen 1994). The interpretation of 
these observations in terms of bulk plasma outflow is not conclusive, however, 
as this effect could be related to the pattern speed of the emitting regions 
rather than to bulk plasma ejection. Another argument for relativistic plasma 
outflow yields a lower limit to the Doppler factor V from the measured radio 
flux density, the angular diameter of the radio emission region, and the upper 
limit of the self-Compton X-ray flux (e.g., Marscher 1987; Ghisellini 1989). An 
accurate measurement of T> through this method requires contemporaneous 
radio and X-ray measurements and, moreover, an accurate determination of 
the radio self-absorption frequency. These conditions are met only rarely, but 
do point to relativistic motions in some flat-spectrum radio quasars. 

Other tests for beaming try to establish conditions for the impossibility of 
intense, rapidly variable emission from stationary radiation sources. The ar- 
gument of Elliot & Shapiro (1974) contrasts the range of allowed black hole 
masses for luminosities governed by Eddington-limited accretion, and variabil- 
ity time scales constrained by the light-travel time across a region with dimen- 
sions corresponding to the Schwarzschild radius of the black hole. A stationary, 
Eddington-limited emitting region is not possible if L4g/At(days) ^> 1, where 
10 n L n ergs s _1 and At(days) are the observed luminosities and variability 
time scales in days, respectively. Klein-Nishina corrections must be applied 
(Dermer & Gehrels 1995) for observations at 7-ray energies. 

Opaqueness of the emitting region to 7-7 attenuation has also been used to 
argue in favor of beaming (Maraschi et al. 1992). In its simplest form, the 
constraint from the compactness parameter £ = 4irm e c 3 /cr T = 6.6 x 10 29 ergs 
s _1 cm -1 implies that if L 45 /At(days) ^> 1, then beaming is implied. Here 
the luminosity refers to emission near 1 MeV where the pair attenuation cross 
section is largest; consequently this test is most sensitive for observations near 
1 MeV. Otherwise, assumptions about the cospatial origin of gamma rays 
and lower energy radiation must be justified by observations of correlated 
variability, since the cross section of >, 100 MeV photons with each other is 
negligible (Dermer & Gehrels 1995). 

Correlated X-ray and TeV observations of Mrk 421 (Macomb et al. 1995; 
Buckley et al. 1996) and Mrk 501 (Catanese et al. 1997) have demonstrated 
that 2-10 keV X-rays and >, 300 GeV 7 rays originate from the same region, 
verifying the cospatial assumption for these sources. Their luminosities are not 
large enough to establish beaming through 7-7 transparency arguments, but 
can be used to determine the mean magnetic field H in the emitting region 
and establish a lower limit to the Doppler factor V through a newly proposed 



2 



beaming test. This test is described in more detail below. 

Because 7-ray observations probe the region nearest the black hole, it is im- 
portant to critically examine these tests. An early hope was that such mea- 
surements could discriminate between accelerating and decelerating jet models 
(e.g., Marscher 1999) by charting the variation of T>, thereby revealing whether 
the evolution of a blazar flare is accompanied by a prompt phase of Doppler 
variation. The possibility that the Doppler factor of the emitting region can 
change, however, introduces intrinsic variability which must be distinguished 
from variability produced by radiative cooling of the emitting particles. It is 
therefore important to consider processes which change the bulk Lorentz fac- 
tor of the radiating plasma. This is done in Section 2. In Section 3, we present 
numerical simulation results showing the effects of bulk plasma (or plasmoid) 
deceleration. Implications for beaming tests and blazar models are discussed 
in Section 4. 



2 Blast- Wave Physics 

Crucial for understanding variability behavior in blazars is to treat the injec- 
tion of relativistic nonthermal particles in the comoving plasma fluid frame 
properly. The blast- wave physics developed to model GRB afterglows (e.g., 
Vietri 1997; Waxman 1997; Wijers et al. 1997) offers a solution to the particle 
injection problem, and provides a method to deal with the deceleration of the 
emitting plasma. The basic idea is that the energy of the injected nonthermal 
particles comes at the expense of the directed bulk kinetic energy of the fluid. 
The variation of the bulk Lorentz factor T of the radiating fluid can be sim- 
ply obtained in a one-zone approximation through a momentum conservation 
equation (Dermer & Chiang 1998). 

Suppose that the system produces an outflow with total energy E and initial 
bulk Lorentz factor r . Because most of the energy of the flow is bound up 
in the kinetic energy of baryons, assumed here to be protons, then E = 
r N th m p c 2 , where N th is the total number of protons. Thus r represents the 
baryon loading of the system. 

It is straightforward to write a conservation equation for the radial (or x) 
momentum component of the fluid, given by 

00 

Tl x (x) = m p cP{(l + a)N th + J dp ■ 7 • [N pi (p; x) + aN e (p; x)]} , (1) 



where a = m e /m p , P = BY = ^/T 2 — 1, and p = (3^ = V7 2 ~~ 1- The func- 
tions N k (p;x) = dN k (p; x) / dp represents the comoving distribution functions 
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of particles of type k = pr (protons) or k = e (electrons) at location x. This 
expression assumes no particle escape. Eq. (1) indicates that the momentum 
of the bulk plasma consists of both the inertia from the thermal protons asso- 
ciated with the baryon loading of the explosion, and the inertia bound up in 
the nonthermal proton and electron distributions. The latter functions evolve 
when nonthermal particles are injected into the plasma and when the energy 
of the nonthermal particles is radiated. 

As a plasmoid or blast wave traverses the surrounding medium, it intercepts 
and sweeps up material. A proton and electron pair is captured by the plasmoid 
with a Lorentz factor V in the comoving plasma frame. The plasmoid captures 
protons and electrons from the surrounding medium at the rate 

dN pi!SW (p,x)/dx = dN e , sw (p, x) / dx = n cxt (x)A(x)5(p - P) , (2) 

The quantity n cxt (x) is the density of particles in the surrounding medium, and 
A(x) is the cross-sectional area of the plasmoid which is effective at sweeping 
up material from the external medium. The power of nonthermal particle 
kinetic energy injected into the comoving frame is simply 

oo 

E ke = m p c 2 J dp (7 — 1) [JVpr )SW (p, x) + aN e>sw (p, x)] , (3) 


where the time derivative refers to time in the comoving frame. The distance 
5x traveled during the comoving time interval St is Sx = St/(BTc). Eqs. (2) 
and (3) therefore imply 

E ke = m p c 2 BT(T - 1)(1 + a)cn cxt {x)A{x) (4) 

(Blandford & McKee 1976). It is important to note that the fraction l/(l+a) = 
99.95% of the energy initially injected into the comoving frame is carried by 
protons, though plasma processes can be effective at transforming this energy 
to electrons or magnetic field. 

By using eq. (1) to write U x (x + Sx), to which is added a term (dTl™ d /dx)Sx to 
account for radiation losses, one obtains an equation of motion for the dynam- 
ics of the blast wave by expanding to first order in Sx and using momentum 
conservation, i.e., IL x (x + Sx) = IL x (x). It is 

-dP(x)/dx = n ext (x)A(x)T(x) 

P(x) (1 + a)N th + JS°dp- 1 - [N pi (p, x) + aN c (p; x)] ' U 

If external Compton scattering processes operate, an additional term must be 
added to take into account the momentum impulse from the scattered external 
photons (Bottcher & Dermer 1999). 

Eq. (5) is the basic equation for calculating the dynamics of a plasmoid by 



4 



sweeping up material from the surrounding medium, and can be solved in a 
number of limiting cases. In the relativistic limit (T ^> 1) and the blast-wave 
case where the area A(x) oc x 2 , there are two important regimes: the adia- 
batic (or non-radiative) and radiative regimes, where the swept-up particles 
retain all or none of their kinetic energy, respectively. Considering the simplest 
case where the density of the external medium can be parameterized by the 
expression n cxt (x) = Uq{x/ 'xd)~ v , one finds that T(x) = Tq for x <C xj, and 
T(x) oc x~ g for x ^> Xd, where g = 3 — i] and g = (3 — i])/2 in the adiabatic 
and radiative regimes, respectively. The deceleration radius 

Xd = ^hnovlmlc^ 3 (6) 

(e.g., Rees & Meszaros 1992), and represents the characteristic distance be- 
yond which the behavior changes from a coasting solution to a decelerating 
solution. The term /& represents the fraction of the full sky into which the 
explosion energy is ejected. 

For intermediate radiative regimes where a fraction ( of the swept-up energy is 
retained in the blast wave, so that a fraction (1 — () is dissipated as radiation 
or lost through escape, Dermer et al. (1999) show that g — (3 — rj)/(l + (). 
In general, one may describe the dynamics of a blast-wave which decelerates 
by sweeping up material from a surrounding medium which is distributed 
according to the relation n cxt (x) oc x~ v by an expression of the form 

r(x) = r /[i + (x/x d y]. (7) 

The relationship between the location x of the blast wave and the observing 
time t bs can be obtained by noting that the radiating element travels a dis- 
tance 5x = VFBcStobs/ (1 + z) during the observing time interval <5t bs, where 
z is the redshift of the source, V = [r(l — £?cos#)] -1 is the Doppler factor, 
and 6 is the angle between the direction of motion of the radiating element 
(or the jet axis) and the observer's direction. From this relation, one can show 
that x oc t Q bs when x <C x,i, and x oc t^ 9+1 ^ when x ^> x^. 

The above expressions are sufficient to treat analytically the basic effects from 
blast wave deceleration and energization through the sweep-up process. As- 
suming that the fraction (1 — () of the swept-up power is dissipated in the 
form of radiation, then the radiated power in the comoving frame at x is 
E oc (1 — ()T 2 (x)n cxt (x)A(x). The received power from a portion of the blast 
wave directed along the line-of-sight to the observer is equal to E amplified 
by a factor T 2 (x)/(1 + z) 2 due to the transformations of energy and time. If 
the observer is outside the Doppler cone of the plasmoid, the emission is weak 
until the radiating plasma has slowed down sufficiently so that the Doppler 
cone intercepts the line of sight. At later times, the emission approaches the 
behavior found in the case where the observer's line-of-sight intercepts the 
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radiating region. 



If ip denotes the opening angle of the plasmoid (or jet), then two limits are im- 
portant when the area of the plasmoid increases oc x 2 . From the above discus- 
sion, the received bolometric power P(t bs) oc (1 — ()T 4 (x)n cxt (x) A(x) / (1 + z) 2 . 
The dynamics of the blast wave changes from a coasting solution to a deceler- 
ating solution when it passes the deceleration radius Xd, which occurs at the 
observing time 

t d = T (l-B co S 9)(l+z)x d /(cT ) -> (1 + 3)^/(21*0) ; (8) 



here B = y 1 — T 2 and the right-hand expression of eq. (8) refers to the case 
when the plasmoid is directed along the line-of-sight. 

For the case 6<ip, P p (t ohs ) oc t 2 ^ 1 when t obs < t d , and P p {t ohs ) oc t^ v ~ 49)/(29+1) 
when t bs 3> t d - For observations at 6 >, ijj, the emission from the blast wave 
begins to intercept the observer's line-of-sight when 6 = l/F(x), which occurs 
when t obs ks t d (r 9)( 29+1 ^ 9 . At this time, the received power is nearly equal to 
the value supposing that 9 < ip. 

We therefore see that the simplest model employing blast-wave physics shows 
how a plasmoid is energized by sweeping up material and converting it into 
nonthermal particle energy in the comoving frame (see Dermer et al. 1999 
for more details). For a jet directed along the observer's line-of-sight (which 
is generally thought of as the standard model for blazars), the sweeping-up 
process produces a flare with bolometric flux rising oc t 2 ^ '. After a sufficient 
amount of material has been swept up to cause the blast wave or plasmoid to 
decelerate, the Doppler deboosting overpowers the additional energization to 
cause the received flux to decay oc t^J'" 49 ^^ 29 ^. For adiabatic (g = 3/2) and 
radiative (g = 3) blast waves in a uniform surrounding medium with r] = 0, 
the light curves decay oc t~£ s and oc t~^^ 7 , respectively. The decaying flux is 
entirely a consequence of the decreasing Doppler factor. Thus it is not valid 
to interpret a decaying flux as evidence for cooling of the emitting particles 
without a discriminant between the effects of cooling and deceleration. 



3 Numerical Calculations 



Using the code described in Chiang & Dermer (1999), we can illustrate the 
effects described in the previous section for a pure synchrotron flare. In the 
simulation shown in Fig. 1, we let the central engine eject 10 48 ergs of plasma 
with a baryon loading given by T = 40 into 10% of the full sky. Thus the 
opening angle of the jet is 11.5°. The jet plasma passes through a medium with 
a uniform density of 0.01 cm -3 , and as it sweeps up this material it converts 
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Fig. 1. Calculations of synchrotron spectra at times given by the labels on the figures. 
Observer is viewing down the jet axis. Parameters of the calculation are given in the 
text. Dashed curve shows the spectrum averaged over a one-day integration time. 

it with high efficiency into nonthermal power-law electrons with a number 
injection index dN/d^i oc 7~ 3 . The energy of the injected electrons ranges 
between 7 = V and 1% of the maximum energy given by balancing the electron 
synchrotron loss time scale and the Larmor time scale in a magnetic field. We 
assume that the magnetic field energy density is 10% of the downstream energy 
density of the swept-up particles. This represents a magnetic field of 0.5 Gauss 
during the phase before the blast wave begins strongly to decelerate. 

Electron synchrotron cooling is taken into account in the calculation, but 
makes only a small contribution to the variability shown here, which is due 
overwhelmingly to energization of the plasmoid by sweeping up particles, and 
to the subsequent deceleration and Doppler deboosting that results from this 
process. The general progress of the flare for the given parameters is to rise 
rapidly at X-ray and soft 7-ray energies. The flare then sweeps to lower energies 
on a much longer time scale. Note that the flare reaches larger vF v peak fluxes 
at higher energies where it is most variable. At lower photon energies, the 
variability is less extreme, and the peak vF v value reached is lower. 

The dashed curve in Fig. 1 shows a one-day time integration over the flaring 
emission. As can be seen, the time-integrated spectrum is much softer than 
the time-resolved spectra because it represents a superposition of hard spectra 
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Fig. 2. Calculations of synchrotron flare light curves using the same parameters 
as in Fig. f. Solid, dashed, and dotted curves are 1 keV X-ray, 1 eV optical, and 
2.4 GHz radio light curves, respectively. Thick curves correspond to observations 
down the axis of the jet with opening angle ip = 11.5°, and thin curves represent 
observations at 6 = 20° from the jet axis. 

which peak at successively lower energies. When fitting blazar flare data using 
a model such as described here, it is therefore necessary to perform time- 
integrations over the flaring emission appropriate to the sampling time of 
the detector. Because gamma-ray telescopes require long observing periods to 
accumulate sufficient statistics, this is especially important for jointly fitting 
hard X-ray/soft 7-ray data, or MeV-GeV flares resulting from the synchrotron 
self-Compton (SSC) or external Compton scattering process. 

Fig. 2 shows light curves for the model synchrotron flare shown in Fig. 1 at X- 
ray, optical and radio frequencies, both along the jet axis (8 = 0°; thick curves), 
and at 20° to the jet axis (thin curves). When observing along the jet axis, the 
peak vF v flux measured at higher frequencies is much greater than the peak 
vF v flux measured at lower frequencies. By contrast, when observing outside 
the opening angle of the jet, i.e., when 9 > ip, one sees that the range of peak 
fluxes becomes much less. Consequently, a flux-limited telescope observing at 
higher photon energies will be much more likely to detect beamed sources 
along the jet axis than a flux-limited telescope sensitive to lower energies, 
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Fig. 3. Total energy received at different observing energies for the flare parameters 
shown in Fig. 1, for observations along the jet axis (solid curve) and at 20° off the 
jet axis (dashed curve). Flux-limited telescopes will preferentially detect jet sources 
along the jet axis when observing at the highest photon energies. 

which will detect on-axis and off-axis sources with nearly equal likelihood. 

This effect is a consequence of the energization and deceleration of the radi- 
ating plasma which causes the received emission to sweep to lower energies at 
late times, and has important implications for any statistical analyses of jet 
sources. In usual statistical treatments, one generally assumes that the relative 
flux observed at different angles to the jet axis is governed by the factor D 3+a , 
independent of photon energy (here a is the energy index of the flux density 
F u oc z/~ Q ; see, e.g., Urry & Padovani 1995). As shown here, the situation is 
much more complicated for flaring sources, which blazars most certainly are. 
A flux-limited X-ray survey primarily detects X-ray jet emission from aligned 
sources, with off-axis jet sources being too faint to detect at X-ray energies. 
By contrast, a flux-limited survey at radio energies will detect on-axis and off- 
axis sources at comparable flux levels. Thus we expect and do see the parent 
population of radio quasars, namely the radio galaxies. The parent population 
of X-ray selected BL Lac objects, by contrast, are at such a low flux level to 
hardly be detectable. 

Fig. 3 shows the integrated energy fluence measured from a blazar flare at 
different photon energies, illustrating the effect just described. We have con- 
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sidered only the synchrotron emission up to this point, but the same behavior 
operates for the SSC flux (Dermer, Mitman, & Chiang 1999, in preparation). 
Variability will be greatest at higher photon energies, where the largest vF v 
fluxes are reached for on-axis sources. The variability time scale will be longer 
and the flux roughly equal for aligned and misdirected blazars at lower photon 
energies. This same general behavior probably applies to external Compton 
scattering emission as well, though additional complications from the narrower 
beaming cone of ECS compared to the synchrotron and the SSC processes 
(Dermer 1995) must be taken into account in this case. 



4 Beaming Tests and Blazar Models 



As noted in the Introduction, the correlated X-ray and TeV observations pro- 
vide a new test for beaming in blazars. This test has been recently applied to 
observations of Mrk 501 by Catanese et al. (1997), Dermer (1998), and Kataoka 
et al. (1999). In this test, it is assumed that the X-ray emission is nonthermal 
synchrotron emission, and that the variability is a result of synchrotron losses 
in a magnetic field of mean intensity H . If the minimum variability time scale 
measured at energy E is denoted by 5t^, then 

(e.g., Tashiro et al. 1995; Takahashi et al. 1996). 

An upper limit on the mean magnetic field H in the emitting region is implied 
because the electrons producing the highest energy synchrotron emission have 
Lorentz factors >2x 10 6 i?c(TeV)(l + z)/T>, where Ec(TeV) is the measured 
energy in TeV of the highest-energy gamma-rays. Synchrotron emission cor- 
related with the TeV flux requires that the electrons radiate in a magnetic 
field at least as great as H(G) = lle ohs>syQ V / [EQ(TeV)(l + z)], where e bs,s y n is 
the measured dimensionless energy of the highest energy synchrotron photons 
produced by the electrons which produce the TeV radiation. When compared 
with the value of H inferred through equations (9), one obtains an expression 
for the Doppler factor, given by 



:i + ^)[E c (TeV)] 3 / 2 
^:(keV)(5tr) 1/2 (^o bs , sy n) 3 / 4 



(Dermer 1998). A lower limit to V is obtained if the TeV flux does not exhibit a 
clear cutoff due to the high-energy cutoff in the electron distribution function. 

With the advent of air Cherenkov telescopes such as Whipple and HEGRA 
detecting BL Lac objects to TeV and tens of TeV energies, this test could in 
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principle provide the largest inferred Doppler factors of all known beaming 
tests (J. H. Buckley, 1998, private communication). It is necessary, however, 
to discriminate clearly between bulk deceleration effects and radiative cool- 
ing effects. Unfortunately, the effects of deceleration mimic those of radiative 
cooling in a variety of ways (Chiang 1999). This is true both for the energy- 
dependent time lags produced by synchrotron cooling, and for the clockwise 
loop diagrams produced by flaring sources when data is plotted in a spec- 
tral index/intensity display (see also Kirk et al. 1999). One such discriminant 
might be the slower decay of the SSC emission compared to synchrotron emis- 
sion (Dermer 1998), but a detailed numerical simulation will be required to 
fully resolve this question. 

Finally, we note that the existence of the process of plasmoid deceleration 
weakens arguments (Buckley 1998) against hadronic models based on the 
long cooling time scales of protons through photo-meson, photon-pion, and 
proton synchrotron processes (e.g., Biermann & Strittmatter 1987) compared 
to the observed rapid variability time scales of BL Lacs. As shown here, nei- 
ther extremely high energy protons nor intense radiation fields are required 
to produce rapid variability, which can result solely from Doppler deboosting 
due to the deceleration of the radiating region. This, and the fact that most 
of the nonthermal particle energy injected into the plasmoid is initially in the 
form of relativistic protons according to the blast wave physics described here, 
seems to give new life to hadronic models of blazars (e.g., Mannheim 1993; 
Gaisser et al. 1995) provided, of course, that hadronic models can successfully 
fit the generic two component vF v blazar spectrum. If radiative cooling does 
not produce the variability behavior, hadronic models must also contend with 
low radiative efficiencies in an uncooled model (see Bottcher & Dermer 1998; 
Totani 1998 for related questions in GRB models). The variability behav- 
ior normally attributed to radiative cooling processes would, in the hadronic 
models, instead be a consequence of plasmoid deceleration. Thus a discrimi- 
nant between cooling and deceleration effects is also of central importance to 
distinguish leptonic and hadronic models of blazars. 



5 Summary 

Recent observations of power-law X-ray afterglow observations of GRBs have 
stimulated the development of new physics for calculating radiation from rel- 
ativistic plasma outflows. The energization of nonthermal particles in the ra- 
diating plasma occurs through a process of sweeping up material from the 
surrounding medium. The internal nonthermal particle energy is extracted 
from the directed kinetic energy of the bulk plasma, causing the bulk plasma 
to decelerate. We have examined idealized flaring behaviors produced when 
bulk plasma sweeps up particles from a uniform medium. Inhomogeneities in 



11 



the surrounding medium as well as in the relativistic plasmoid will complicate 
the situation. 

We have shown that 

• Bulk plasma deceleration effects must be included in models and statistical 
studies of blazars; 

• Time integrations over the varying Doppler factor of the radiating plasma 
must be performed when fitting blazar data, employing integration ranges 
appropriate to the observing times of the different telescopes; 

• A newly developed test for beaming in blazars using correlated X-ray and 
TeV variability observations must answer the criticism that the variability 
is due not to radiative cooling but rather to Doppler deboosting; 

• Hadronic models for blazars do not have to confront the difficulty of the long 
radiative cooling time scales of hadrons, since variability can be achieved 
through plasmoid deceleration. 

Important for further progress on these questions is an observational discrim- 
inant between cooling and deceleration processes. 
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